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In ultra-central heavy-ion collisions, anisotropic hydrodynamic flow is generated by density fluc¬ 
tuations in the initial state rather than by geometric overlap effects. For a given centrality class, the 
initial fluctuation spectrum is sensitive to the method chosen for binning the events into centrality 
classes. We show that sorting events by total initial entropy or by total final multiplicity yields event 
classes with equivalent statistical fluctuation properties, in spite of viscous entropy production dur¬ 
ing the fireball evolution. With this initial entropy-based centrality definition we generate several 
classes of ultra-central Pb+Pb collisions at LHC energies and evolve the events using viscous hy¬ 
drodynamics with non-zero shear but vanishing bulk viscosity. Comparing the predicted anisotropic 
flow coefficients for charged hadrons with CMS data we find that both the Monte Carlo Glauber 
(MC-Glb) and Monte Carlo Kharzeev-Levin-Nardi (MC-KLN) models produce initial fluctuation 
spectra that are incompatible with the measured final anisotropic flow power spectrum, for any 
choice of the specific shear viscosity. In spite of this failure, we show that the hydrodynamic model 
can qualitatively explain, in terms of event-by-event fluctuations of the anisotropic flow coefficients 
and flow angles, the breaking of flow factorization for elliptic, triangular and quadrangular flow mea¬ 
sured by the CMS experiment. For elliptic flow, this factorization breaking is large in ultra-central 
collisions. We conclude that the bulk of the experimentally observed flow factorization breaking 
effects are qualitatively explained by hydrodynamic evolution of initial-state fluctuations, but that 
their quantitative description requires a better understanding of the initial fluctuation spectrum. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 


I. INTRODUCTION 

Heavy-ion collisions at the Relativistic Heavy-Ion Col¬ 
lider (RHIC) and the Large Hadron Collider (LHC) of¬ 
fer a unique window to study the physics of strongly 
interacting matter in the quark-gluon plasma (QGP) 
phase. An interesting fundamental property of the QGP 
is its viscosity which, in combination with the expan¬ 
sion rate of the fireball created in the heavy-ion colli¬ 
sion, controls its fluidity. Specifically, the shear viscos¬ 
ity to entropy density ratio r\/s controls the efficiency 
with which spatial inhomogeneities and anisotropies of 
the initial pressure gradients in the fireball are converted 
to anisotropies in the final collective flow pattern |Tj. 
Spatial anisotropies in the initial pressure distribution 
(related to the initial energy density profile through the 
equation of state (EOS)) are typically characterized in 
terms of a set of initial eccentricity coefficients e n (de¬ 
fined below as Fourier coefficients of the initial energy 
density distribution in the plane transverse to the beam 
direction), while the final flow anisotropies are quantified 
through coefficients v n obtained from a Fourier decom¬ 
position of the final transverse momentum distribution 
of the emitted hadrons. Large shear viscosity tends to 
smooth out flow anisotropies and thus suppress the con¬ 
version efficiency v n /En BUS- 

The experimental observation of large flow anisotropies 
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at both RHIC @3 and LHC 5j indicates that the QGP 
has very low specific shear viscosity 77 /s P] close to its 
quantum limit Eld. However, since the extraction of 
r\ls from experimental data goes through the conversion 
efficiency v n /e n [5j, it requires, in addition to a mea¬ 
surement of v n , a knowledge of the initial eccentricity 
spectrum s n . As the initial density profile of the fireball 
created in the collision fluctuates from event to event, so 
do its initial eccentricity and final flow power spectra, e n 
and v n . At finite impact parameter the eccentricity spec¬ 
trum e n receives contributions from both collision geom¬ 
etry and event-by-event density fluctuations. In ultra¬ 
central collisions, on the other hand, overlap geometry 
and associated model uncertainties can be largely elim¬ 
inated as a contributing factor ]9], focussing attention 
on the initial spectrum of quantum fluctuations of the 
strongly interacting quantum fields within the colliding 
nuclear wave functions. These are notoriously difficult 
to compute theoretically, making experimental studies 
of ultra-central heavy-ion collisions hhhh particularly 
valuable for improving our quantitative understanding of 
heavy-ion collision dynamics. 

It bears pointing out that a measurement of the fi¬ 
nal flow power spectrum v n in ultra-central collisions be¬ 
tween equal-size nuclei provides access to both this initial 
fluctuation spectrum [l.'fl iHj (the “Little Bang temper¬ 
ature power spectrum” [15]) and the QGP transport co¬ 
efficients, especially its shear viscosity mug ns]. Since 
the specific shear viscosity 77 /s controls the anisotropic 
flow response of the expanding liquid to the initial den¬ 
sity fluctuations, and neither the specific shear viscos¬ 
ity nor the initial fluctuation spectrum are theoretically 
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FIG. 1: (Color online) The root mean square (rms) initial eccentricity e n {2} for n = 2—5, as a function of centrality, for different 
centrality selection criteria: using the number of participant nucleons 7V part (solid black), the number of binary collisions lV co n 
(red dashed), impact parameter b (green dotted), or initial total entropy dS/dy (dotted blue). For this figure 1 million events 
were generated from the MC-Glauber model with minimum inter-nucleon distance (“hard core radius”) r m j n = 0.9fm [191 }20] . 
for Pb+Pb collisions at = 2.76 A TeV, with binary collision to wounded nucleon ratio a = 0.118. 


known with any precision, a phenomenological analysis 
of anisotropic flow fluctuations must be sufficiently com¬ 
prehensive to be able to constrain both simultaneously. 

The CMS Collaborations at the LHC has measured 
023 the anisotropic flow coefficients of charged hadrons 
in ultra-central Pb+Pb collisions at \/s = 2.76 A TeV, us¬ 
ing events from the tail of the charged multiplicity distri¬ 
bution corresponding to 0—0.2% centrality. Implement¬ 
ing such tight multiplicity cuts is a challenge for the¬ 
oretical simulations which have typically much smaller 
event samples than the experimental measurements. In 
Sec. [Tl] of this work, we investigate how different ways of 
cutting centrality bins can affect the initial-state eccen¬ 
tricity spectrum, and show how to faithfully mimic the 
experimental centrality selection procedure without the 
need for evolving astronomically large numbers of initial 
configurations. In Sec. |III[ we show the charged hadron 
particle spectrum and its anisotropic flow coefficients. In 
Sec. m we check the flow factorization ratios in ultra¬ 
central Pb+Pb collisions at the LHC. We summarize our 
findings in Sec. [V] 

II. INITIAL STATE FLUCTUATIONS AND 
EVENT SELECTION IN ULTRA-CENTRAL 
NUCLEUS-NUCLEUS COLLISIONS 

Experimentally, the ultra-central collision events are 
defined as those with the highest measured final charged 


multiplicities. However, in theoretical calculations, using 
final charged multiplicities to determine the event cen¬ 
trality is numerically very expensive. Especially for the 
narrow 0-0.2% centrality bin, the numerical effort spent 
on evolving the remaining 99.8% of the simulated events 
is wasted. So we would like to determine (at least es¬ 
timate) the event centrality of each initial profile before 
evolving them through viscous hydrodynamics. 

We begin by showing that for ultra-central collisions 
the initial fluctuation spectrum {e n {2}\n = 1, 2 ,... } 
(whose precise definition is given further below) can de¬ 
pend strongly on the method for determining centrality. 
We study four popular choices of the collision parame¬ 
ters (fVp ar t, A co u, 6, dS/dy) to categorize the event cen¬ 
trality of every initial density profile generated from the 
MC-Glauber model [22]. The produced initial entropy 
density at To in the transverse plane is calculated as 

s(rj_;To) = ^nwN(u) +a»Bc(i’i)j , (1) 

where n is an overall normalization factor, tuned to re¬ 
produce the final charged hadron multiplicity dN ch /dii 
at mid-rapidity, and a is the mixing ratio between the 
wounded nucleon (nwN( r )) and binary collision (nBc( r )) 
profiles [22]. We choose a = 0.118 [25] to repro¬ 
duce the measured centrality dependence of dN ch /dri in 
2.76 A TeV Pb+Pb collisions. 

In Fig. [l] the initial eccentricities s n {2} (n = 2—5) are 
shown as functions of centrality. For semi-peripheral and 
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( 0~0.2%) dN / d y 

(0-0.2%) ds/dy 

(0~l%)dN/dy 

(0-l%) d s/dy 

(1-2 %) d N/dy 

(1-2 %)dS/dy 

( 2-3%) d N/dy 

(2-3 %) dS /dy 

£2{2} 

0.082 + 0.002 

0.081 + 0.002 

0.085 + 0.001 

0.086 + 0.001 

0.096 + 0.001 

0.096 + 0.001 

0.114 + 0.001 

0.114 + 0.001 

es{2} 

0.084 + 0.002 

0.085 ± 0.002 

0.088 + 0.001 

0.088 + 0.001 

0.092 + 0.001 

0.092 + 0.001 

0.097 + 0.001 

0.096 + 0.001 

e 4 {2} 

0.095 + 0.002 

0.096 + 0.002 

0.095 + 0.001 

0.095 + 0.001 

0.097+ 0.001 

0.097 + 0.001 

0.102 + 0.001 

0.103 + 0.001 

£5(2} 

0.089 + 0.002 

0.092 + 0.002 

0.097 + 0.001 

0.097 + 0.001 

0.104 + 0.001 

0.104 + 0.001 

0.110 + 0.001 

0.110 + 0.001 


TABLE I: Comparison of the initial eccentricities e n {2} (n = 2—5), computed with the MC-KLN model with a minimum 
distance (“hard core”) r m i n = 0.4fm between nucleons | 211I41| . in central and ultra-central centrality bins that are determined 
by initial total entropy or final charged hadron multiplicity, respectively. 


peripheral collisions, the e„{2} spectrum is insensitive to 
the method used to determine the centrality bins. How¬ 
ever, differences become noticeable when we select events 
within the top 1% centrality range. We notice that using 
the initial total entropy to determine the centrality bin 
gives the largest initial eccentricities {e„{2}}, about 10% 
larger than if we use the number of participant nucleons. 
So for a narrow bin of very central collisions, extra care 
is needed in the event selection. 

Experience based on the hydrodynamic evolution of 
smooth, ensemble-averaged initial conditions has shown 
that the observed charged multiplicity in the final state, 
dNch/drj, used by the experiments for cutting on central¬ 
ity, is monotonically related to the initial entropy density 
dS/dy. Figure [2] demonstrates that this remains true on 
average, within a narrow variance caused by fluctuations 
of the viscous entropy production in fireballs with par¬ 


ticularly rough initial conditions (examples are shown in 
Fig.§, when one evolves bumpy initial conditions that 
fluctuate from event to event. 1 Note that the events plot¬ 
ted in Fig.[2]are all within the 0—5% centrality bin where 
variations in dS/dy are caused almost entirely by fluctu¬ 
ations, and variations in overlap geometry play only a 
small role. 

Table |T] shows that, within the statistical uncertainties 
imposed by our limited number of fully evolved events, 
the eccentricity coefficients £ 2 , 3 , 4,5 defined in Eq. ([ 2 ]) be¬ 
low for events in a given centrality bin agree for bins de¬ 
fined by ordering the events according to initial dS/dy or 
final dN c h/dy. We also checked and confirmed that the 
distributions of these coefficients within the bins agree 
for both centrality definitions. Based on all these ob¬ 
servations, we proceeded to generate 1 million minimum 
bias Pb+Pb events as described next, ordered them by 
their initial total entropy dS/dy , divided them into cen¬ 
trality classes (as shown in Fig. [ 2 ]) and then took the 
0.2% of them with the largest dS/dy values to represent 
the 0-0.2% centrality class of ultra-central Pb+Pb colli¬ 
sions studied by the CMS Collaboration. These events 
are then used in the rest of the paper. 

In order to capture the bias in ultra-central collisions 
towards events whose entropy production fluctuates up¬ 
ward from the mean, we implement pp multiplicity fluc¬ 
tuations in the MC-Glauber model in a way that allows 
us to reproduce the experimentally observed KNO scal¬ 
ing of the pp multiplicity distributions. Our implemen¬ 
tation is described in detail in [22] (see earlier work in 
f27M29| for related but slightly different approaches). 2 
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FIG. 2: (Color online) Correlation between the initial to¬ 
tal entropy dS/dy and final charged multiplicity dN c ^/dy for 
(ultra-)central Pb+Pb collisions at the LHC (0-5% central¬ 
ity). The plot is based on 10,000 previously generated viscous 
hydrodynamic events with fluctuating MC-KLN initial con¬ 
ditions using a minimum inter-nucleon distance r m i n = 0.4 fm 
mum and evolved with 77 /s = 0.2. The lines separate differ¬ 
ent centrality classes (as indicated), defined by ordering the 
events according to dN c ^/dy (horizontal lines) or dS/dy (ver¬ 
tical lines). The normalized variance of the distribution is 
about 0.6% in both vertical and horizontal directions. 


1 Fig- I is based on events with MC-KLN initial conditions that 
were evolved with rj/s = 0 . 20 , the largest value for the specific 
shear viscosity studied in this work. For MC-Glauber initial con¬ 
ditions evolved with 77 /s = 0.08 the viscous entropy production 
is smaller, and we correspondingly found a reduced normalized 
variance of 0.2% for the analogous dS/dy vs. dN^/dy distribu¬ 
tion. Both sets of hydrodynamic events were previously gener¬ 
ated for an earlier study EE| and do not include pp-multiplicity 
fluctuations which, as discussed below, tend to increase the 
“bumpiness” of the initial density profiles and are therefore ex¬ 
pected to somewhat increase the normalized variance of the dis¬ 
tribution shown in Fig.|2l 

2 For the MC-KLN model we currently do not implement any 
collision-by-collision multiplicity fluctuations. An implementa- 
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FIG. 4: (Color online) The initial total entropy distribution 
(a) and the initial rms eccentricity £n{2} as a function of har¬ 
monic order n (b) from the MC-Glauber and MC-KLN mod¬ 
els for 0-0.2% ultra-central Pb+Pb collisions at 2.76 +TeV. 
The results shown with open symbols are calculated requir¬ 
ing a minimum inter-nucleon distance (“hard core radius”) 
v min = 0.9 fm 1 19] 1201 [22] . The filled symbols in panel (b) are 
obtained by sampling nucleon configurations that incorporate 
a realistic repulsive 2-body nucleon-nucleon correlation [23 1 . 


FIG. 3: (Color online) Initial entropy density profiles in the 
transverse plane from the MC-Glauber model, for one selected 
sampling of the nucleon positions inside the colliding nuclei. 
Panel (a) does not contain pp multiplicity fluctuations. Pan¬ 
els (b) and (c) correspond to two different samplings of the 
entropy production associated with each nucleon-nucleon col¬ 
lision when pp multiplicity fluctuations are included as de¬ 
scribed in the text. 


Figure [3] illustrates the effect of pp multiplicity fluctu¬ 
ations on the shape of initial entropy density profile in 
Pb+Pb collisions. They tend to increase the variance 


tion in the MC-KLN model of KNO scaling in pp coolisions can 
be found in Ref. ESI- 


of the density fluctuations in the initial state, amplify¬ 
ing the magnitude of hot spots and deepening the val¬ 
leys between them. Overall the initial density distribu¬ 
tions become more bumpy without, however, changing 
the characteristic radius of the hot spots (which is still 
given by the nucleon size). In Fig. [4^, we show the ini¬ 
tial total entropy distribution for MC-Glauber and MC¬ 
KLN initial conditions, for Pb+Pb collisions in the 0- 
0.2% centrality bin at LHC energies. Even in the absence 
of pp multiplicity fluctuations, the MC-Glauber model 
yields a significantly broader multiplicity distribution in 
Pb+Pb collisions than the MC-KLN model. By adding 
pp multiplicity fluctuations to the MC-Glauber model, 
the chances for upward fluctuations in the multiplicity 
for ultra-central Pb+Pb collisions are strongly increased. 
For a correct simulation of ultra-central Pb+Pb initial 
conditions, proper inclusion of pp multiplicity fluctua¬ 
tions is therefore mandatory. 
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The shape fluctuations of the initial conditions can be 
characterized by a series of complex eccentricity coeffi¬ 
cients. For a single event and for n > 2, the n-th order 
eccentricity of the initial condition is defined as the mod¬ 
ulus of the complex quantity 

_ in<s> n __ fd 2 rr n e in ^e(r, </>) 
n f cPrr n e(r, <fi) ’ 

in which we use the transverse local rest frame energy 
density e(r, </>) as a weight function. The phase <&„ of the 
complex quantity in Eq. § defines the n-th order par¬ 
ticipant plane angle. Since the energy density is centered 
at the origin, Eq. (|2j) gives zero for n = 1. Hence, one 
defines using the next order in the cumulant expansion 

m, 


Model 

dN/dp Ii^ko.s 

(p T ) (GeV) 

MC-Glb. ideal 

1793.85 ± 1.02 

0.715 ± 0.001 

MC-Glb. p/s = 0.08 

1799.96 ± 1.07 

0.705 ± 0.001 

MC-Glb. p/s = 0.12 

1780.44 ± 1.03 

0.704 ± 0.001 

MC-Glb. p/s = 0.20 

1797.76 ± 1.68 

0.711 ± 0.001 

MC-KLN ideal 

1809.78 ± 0.60 

0.688 ± 0.001 

MC-KLN p/s = 0.08 

1807.66 ± 0.43 

0.682 ± 0.001 

MC-KLN p/s = 0.12 

1807.81 ± 0.43 

0.685 ± 0.001 

MC-KLN p/s = 0.20 

1811.03 ± 0.53 

0.692 ± 0.001 


TABLE II: The total yield of charged hadrons and their mean 
Pt in 0-0.2% in Pb+Pb collisions at the LHC. In each set 
of runs we fixed the overall normalization factor to the ex¬ 
perimentally measured charged hadron multiplicity at 0-5% 
centrality. 


Jd 2 rr 3 e^e(r, <j>) 
f d?r r 3 e(r, </>) 


III. PARTICLE SPECTRA AND THEIR FLOW 
(3) ANISOTROPIES 


s n {2} = is the root mean square value of these 

eccentricity coefficients for an ensemble of fluctuating 
events (where (...) denotes an ensemble average). 

In Fig. we plot the initial fluctuation power spec¬ 
trum e„{2} for the MC-Glauber and MC-KLN models in 
0-0.2% ultra-central collisions at the LHC. At this cen¬ 
trality, the e n {2} are entirely dominated by fluctuations 
in the initial density profile, with no contributions from 
geometric overlap effects. With the r n weighting factor 
in Eq. ([2]), the magnitudes of e„{2} are all similar (ex¬ 
cept for £i{ 2} which is ~ 40% smaller), for both initial 
condition models. Without pp fluctuations, the e n {2} 
from the MC-KLN model are close to those from the 
MC-Glauber model for n > 3 and about 10-30% smaller 
than the MC-Glauber eccentricities for n = 1 and 2. In¬ 
clusion of pp fluctuations in the MC-Glauber model leads 
to a significant increase of e n {2} for all harmonics, rang¬ 
ing from 15% for n = 1 to almost 50% for n = 9. We also 
include results (shown by filled symbols) which sample 
realistic nucleon configurations of a lead nucleus in the 
presence of repulsive two-nucleon correlations |23j . For 
both the MC-Glauber and MC-KLN models, realistic 2- 
body nucleon-nucleon correlations result in low-order ec¬ 
centricities that are very close to those obtained when 
simply imposing (as suggested in m uni) a minimum 
inter-nucleon distance (“hard core radius”) r m i n = 0.9fm 
in the Monte Carlo sampling of the nucleon positions 
from nuclear density distribution. Higher-order e n with 
n > 6 are about 10% larger for samples containing realis¬ 
tic nucleon-nucleon correlations. 

In the next section, we will evolve these initial con¬ 
ditions using viscous hydrodynamics with four different 
specific shear viscosities, p/s = 0,0.08,0.12,0.20. For 
each set of runs we generate 1000 events which is suffi¬ 
cient for a statistical analysis of anisotropic flows. 


In Fig. [5] we show the transverse momentum spectra of 
charged hadrons from the MC-Glauber and the MC-KLN 
models, hydrodynamically evolved with different p/s val¬ 
ues. In ultra-central collisions, we find that the shear vis¬ 
cosity has only minor effects on the slope of the charged 
hadron spectra. 



FIG. 5: (Color online) Charged hadron pr-spectra from the 
MC-Glb and MC-KLN models, for different specific shear vis¬ 
cosities as indicated, for Pb+Pb collisions of 0-0.2% centrality 
at yfs = 2.76 A TeV. 

The MC-Glauber initial conditions, which include 
pp multiplicity fluctuations, result in slightly flatter 
(harder) charged particle spectra compared to the MC- 
KLN model. This is because the additional collision-by¬ 
collision multiplicity fluctuations increase the bumpiness 
of the initial density profiles and thereby the initial pres¬ 
sure gradients, which drives the system to develop more 
radial flow. In Table |TT| we summarize for the eight cases 
shown in Fig. [5] the total charged hadron yields and their 
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FIG. 6: (Color online) pr-integrated v„{2} of charged hadrons in 2.76 A TeV Pb+Pb collisions at 0-0.2% centrality, from 
the MC-Glauber (a,c) and the MC-KLN models (b,d) models. In (a,b) repulsive hard core correlations among nucleons are 
implemented by simply requiring a minimum inter-nucleon distance r m j n = 0.9 fm when sampling the nucleon positions, in (c,d) 
realistic 2-nucleon correlations |23| are included in the sampling routine. The solid squares show the CMS measurements |12| 
for comparison. As in the CMS measurements, v n {2} is integrated over px from 0.3 to 3 GeV. 
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mean px in 0-0.2% ultra-central Pb+Pb collisions at the 
LHC. Ideal hydrodynamic simulations result in a larger 
mean px than the viscous ones. Due to the lack of viscous 
entropy production in an ideal fluid, the system needs to 
start with a higher peak temperature in order to produce 
the same final multiplicity. This results a larger pressure 
gradient and a ~ 10% longer fireball lifetime, which both 
lead to the development of stronger radial hydrodynamic 
flow compared to viscous evolution. Because of the large 
system size in Pb+Pb collisions, viscous effects on the 
fireball expansion are not as pronounced as in smaller 
collision systems. A growth of rj/s by a factor 2.5 is seen 
to increase the charged hadron mean px by just over 1%. 

Due to fluctuations, the px distributions of individual 
events are azimuthally anisotropic even for zero impact 
parameter. These anisotropies can be characterized by 
a series of (px and rapidity dependent) harmonic flow 
coefficients v n and associated flow angles 4+: 

dN _ dN 
dyp T dp T d(j)p dypxdpx 

OO \ 

1 + 2 v n{PT, y) cos (cj)p-^ n (px, y)) J • (4) 

n=l ' 


In the absence of non-flow effects, the “2-particle cumu- 
lant flow” v n {2} is the root mean square of the fluctuat¬ 
ing flow v n in an event ensemble: v n {2} = \J(v%). 

In Fig. [6] we show the charged particle p^-integrated 
v n {2} as a function of its harmonic order n for both the 
MC-Glauber and the MC-KLN models, for 0-0.2% ultra¬ 
central Pb+Pb collisions at the LHC. The MC-Glauber 
initial conditions account for pp multiplicity fluctuations. 
The top and bottom panels compare calculations based 
on two different approaches for including repulsive 2- 
nucleon correlations within the nuclei, as described in 
the figure caption. The differences between them are 
seen to be small and within the experimental uncertain¬ 
ties of the measured py-integrated v n {2} values. For 
the pT-differential flow observables discussed below, the 
corresponding differences are typically less than the line 
widths in the plots. In the rest of the paper we there¬ 
fore show only calculations based on the more realistic 
modeling of repulsive NN-correlations described in [22]. 

Since the initial eccentricities are small on average 
(e„{2}~ 0.07-0.12 for MC-KLN and e n {2} ~ 0.12-0.17 
for the MC-Glauber model), non-linear mode couplings 
[32] in the higher order (n > 3) anisotropic flow coeffi¬ 
cients v n , involving the product of multiple e n [331 :3U, 
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FIG. 7: (Color online) MC-Glauber model calculations, including realistic NN correlations in the initial state, of the pr- 
differential two-particle cumulant v n {2} (n G [1,6]) of charged hadrons, compared with CMS measurements for 2.76 A TeV 
Pb+Pb collisions at 0-0.2% centrality 112] . 


are suppressed. We expect a dominantly linear response 
between the initial e n {2} and the final i+{2}. ! The con¬ 
version efficiency v n {2}/e n {2] is controlled by the spe¬ 
cific shear viscosity of the medium, as seen in Fig. [6] 
Whereas all simulations capture the general trend of in¬ 
creasingly strong suppression of v n for higher harmonics 
n, which is caused by the non-zero shear viscosity of the 
medium, both initial condition models fail to quantita¬ 
tively reproduce the measured v n {2} spectrum (i>„{2} as 
a function of n). This statement holds for all the 77 /s 
values explored in the simulations. 

In the CMS data, the magnitude of v 2 {2} is very 
close to 7; 3 {2}, which can not be reproduced, nei¬ 
ther with MC-Glauber nor with MC-KLN initial con¬ 
ditions. 3 4 In hydrodynamic simulations, the conversion 
efficiency from initial e n {2} to final v n {2} is expected 
to decrease with increasing order n 1 3 GS, roughly as 


3 A recent analysis Eg shows that even in central collisions there is 
a subleading contribution to the triangular flow V 3 that is driven 
by a different radial moment £3 5 of the initial triangularity. Still, 
the response to the linear superposition of these two dominant 
drivers is linear Eg- 

4 The authors of Refs. |291 I36| argue that accounting for short- 
range correlations between nucleons in the Monte Carlo sampling 
of the nucleon positions within the colliding nuclei and for bulk 
viscous effects somewhat alleviates this problem. However, the 
experimentally observed rapid change of the ratio U 2 { 2 }/u 3 { 2 } in 
ultra-central multiplicity bins of increasing width, from a value 
«1.0 at 0-0.2% centrality to «1.15 at 0-2.5% centrality, rising 
to almost 1.5 at 2.5-5% centrality m, remains unexplained. 


\w[v n {2}/e n {2}} oc — (?n 2 [37]. Given this expectation and 
the hydrodynamic results shown in Fig. 9 of Ref. [3] and 
Fig. [6]here, an observed ratio of v 2 {2}/vz{2} ~ 1 in ultra¬ 
central Pb+Pb events would require a ratio of the corre¬ 
sponding initial eccentricities £2{2}/e 3 {2} ~ 0.5—0.7, de¬ 
pending on 77 /s. This is clearly inconsistent with the 
predictions from the MC-Glauber and MC-KLN models 
shown in Fig. [4jr. Starting with an initial fluctuation 
spectrum featuring £ 2 {2}~£ 3 {2}, as shown in Fig. Eh 
for reasonable values of 77 /s the final v 2 {2} from hydrody¬ 
namic simulations will always be about 30% larger than 
t> 3 {2}. Any theoretical calculation using such initial con¬ 
ditions will therefore overestimate the measured v 2 /v 3 
ratio (at least as long as only shear viscous effect are in¬ 
cluded). Within the current hydrodynamic framework, 
to reproduce a v n hierarchy similar to the one observed 
by CMS, we need an initial condition model that pro¬ 
vides a triangular deformation £ 3 {2} that is significantly 
larger than £ 2 {2} in these ultra-central collisions. 

Comparing the overall magnitudes of the predicted 
v n {2} values with the CMS data |12] in Fig. [6] we find 
that the data prefer a relatively large specific shear vis¬ 
cosity, for both MC-Glauber ( 77 /s ~ 0.3) and MC-KLN 
( 77 /s ~ 0.2) initial conditions. In the MC-Glauber case, 
this conclusion differs from that drawn from earlier stud¬ 
ies that did not include pp multiplicity fluctuations. The 
increased initial r.rn.s. eccentricities caused by multiplic¬ 
ity fluctuations generate larger r.rn.s. flow anisotropies 
unless tempered by a larger shear viscosity. Since our 
MC-KLN initial conditions currently do not yet include 
the effects of pp multiplicity fluctuations, the 77 /s values 
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FIG. 8: (Color online) Similar to Fig. [t] but for MC-KLN initial conditions. 


required by that model may further increase once this 
model deficiency is repaired. 

Next, we study the pr-differential anisotropic flows, 
v u (pt)- Since for ultra-central collisions the anisotropic 
flow coefficients are dominated by the fluctuating bumpi¬ 
ness of the initial energy density distributions, the asso¬ 
ciated flow angles are randomly oriented relative to the 
(very small) impact parameter and correlated with each 
other only through the density fluctuations and not by 
overlap geometry. In such fluctuation-dominated situa¬ 
tions, the flow angles T 7) (p T ) feature particularly strong 
PT-dependence, oscillating around their p^-averaged val¬ 
ues in patterns that strongly fluctuate from event to 
event [38]. Different flow observables correspond, in 
general, to different moments of the probability distri¬ 
bution that governs the event-by-event fluctuations of 
v n and The specific flow extraction method used 

in the CMS measurements determines the following pj- 
differential analog of the two-particle cumulant flow of 
charged hadrons: 


w{2 }(pt) 


(Vn{PT)Vn f COs[n(^» (p T ) ~T^ of )]) 

'•;r{2i 


(5) 


Here rF ef {2} and are the px-integrated n-th order 
flows obtained from a set of reference particles, taken 
from a different sub-event (usually a neighboring rapid¬ 
ity bin) to eliminate non-flow and self-correlations. As 
done in the CMS analysis, we choose all charged hadrons 
with transverse momenta between 1 and 3 GeV as the 
reference particles; in the experiment, this choice opti¬ 
mizes the sensitivity for the measurement of higher order 
v n . Since the theoretically calculated momentum distri¬ 
butions are free of non-flow effects and continuous, corre¬ 


sponding to effectively infinite statistics in a single event, 
self correlations can be neglected, avoiding the need for 
different sub-events for signal and reference particles. 

In Figs. [7] and [8] we show the py-differential flows 
v n {2 }(pr) at 0-0.2% centrality from the MC-Glauber and 
MC-KLN models, compared with the CMS data. We find 
that the direct flow v\{2} is insensitive to the p/s value 
used in the simulations. It flips from negative to positive 
at pt ~ 1.2 GeV due to global momentum conservation. 
Comparing zj 2 {2}(pt) through vq{2 }(pr) with the CMS 
data H2], our results using the MC-Glauber model with 
p/s = 0.20 provide a fairly good description, except for 
the differential elliptic flow V 2{2 }{pt) which is overesti¬ 
mated by the calculation. The MC-KLN initial condi¬ 
tions with p/s = 0.20 can describe V 2 { 2 }(pt) better but 
underestimate v 3 {2}(p T ). Although neither model de¬ 
scribes all the data equally well, the overall picture seems 
to favor a QGP shear viscosity in the hydrodynamic sim¬ 
ulations that lies at the upper end of the explored range, 
( v/ s )qgp ~ 0.20. 


IV. FLOW FACTORIZATION 


In a single event, hydrodynamic flow effects on two- 
particle correlations factorize into a product of single¬ 
particle flow coefficients f3£ 40j ■ Due to event-by-event 
flow fluctuations this factorization is broken in mea¬ 
surements that are based on event averaged observables 
mm- The factorization breaking effects can be quan- 
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p“ g -p£ sso (GeV) p^ ig ~Pt s ° (GeV) p“ g -p^ s0 (GeV) p“ g -p£ sso (GeV) 

FIG. 9: (Color online) Flow factorization ratio r n (p-/' s : Pr s °) (n = 2,3,4) from model calculations with MC-Glauber (a-d) and 
MC-KLN (e-h) initial conditions, ultra-central 2.76 71 TeV Pb + Pb collisions at 0-0.2% centrality. Four values of p/s were 
explored in each case, as indicated. Theoretical results are compared with measurements by the CMS collaboration m- 


p“ g e[1.0,1.5] GeV p“ g e[1.5,2.0] GeV p“ g e[2.0,2.5] GeV p“ g e[2.5,3.0] GeV 



p“ g -p^ so (GeV) p“ g -p^ ss0 (GeV) p“ g -p^ s0 (GeV) p“ g -p|? so (GeV) 

FIG. 10: (Color online) Similar to Fig. [9j but for n = 3. 
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FIG. 11: (Color online) Similar to Fig. [9] but for n = 4. 
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tified by the “flow factorization ratios” 

Kia(PT1,PT2) 


r n (pri,PT2) ■ = 


■\/V n A(pTl,PTl)V n A{pT2,PT2) 
(^n(PTl)^n(PT2)cOs[?r(^»(pTl)-^n(PT2))]) 
V( V n(PTl))(vl{pT2)) 


( 6 ) 


where I 4 ia(pti>PT 2 ) are the two-particle differential flow 
coefficients [3Hll40j characterizing the correlation between 
a “trigger particle” at transverse momentum p!/ Ig = pxi 
and an “associated particle” from the same event at 
^ass° _ p T2 _ d ev i a tion of this ratio from one indi¬ 
cates that flow factorization is broken. In fact, if r n is 
larger than one, factorization must be broken by non-flow 
effects - hydrodynamic flow fluctuations always lead to 

r n < 1 EDI- 

In Figs. [9] - [TT] we show the flow factorization ratios 
^ 2 , 3,4 at 0-0.2% centrality, for different values of 77 /s. For 
the given ranges of trigger particle transverse momentum, 
these ratios are shown as a function of p^ lg —p|? so and 
compared with the CMS measurements [T2j . 

Figure [9] reveals that in ultra-central Pb+Pb collisions 
the ratio 7*2 exhibits a large breaking of flow factorization 
at large pp Ig —p|? so . This breaking is not due to the onset 
of non-flow effects at large pr as originally suspected [39], 
but due to flow fluctuations caused by initial-state fluc¬ 
tuations and very well described by hydrodynamic evolu¬ 
tion of the latter. We have checked that about half of the 
factorization breaking arises from the flow angle fluctua¬ 
tions [35] while the other half comes from fluctuations of 
the magnitudes the anisotropic flows. The figure shows 
that 77 /s affects the ratio 7~2 non-monotonically, and we 
found that this is associated with a change in the rela¬ 
tive contribution of v n and fluctuations to the break¬ 
ing of flow factorization. This may indicate a non-trivial 
interplay of viscous damping effects on the fluctuations 
in flow magnitude and flow angle. Indeed, our studies 
showed that increasing the shear viscosity suppresses v n 
fluctuations but strengthens flow-angle correlations [Ill- 

Comparing the upper panels with the lower ones in 
Fig- m we find that, irrespective of the choice of 77 /s, 
the ratio ?~2 deviates from 1 more strongly when we use 
MC-KLN initial conditions than for MC-Glauber profiles. 
Similar statements, but with the opposite sign, hold for 
7~3 and 7 - 4 , shown in Figs. [TO] and |TT] in their case, the fac¬ 
torization breaking effects are weaker for MC-KLN initial 
conditions than for MC-Glauber ones, again with little 
sensitivity to 77 /s. This implies that r n actually responds 
more strongly to changes in the initial fluctuation spec¬ 
trum than to variations of the shear viscosity. The CMS 
i >2 factorization breaking data in ultra-central collisions, 
shown in Fig. [9] favor MC-Glauber-like initial-state fluc¬ 
tuations over the MC-KLN ones, but Figs. [To] and [lT] 
lead to the opposite conclusion: while the flow factoriza¬ 
tion breaking effects for triangular ( 7 - 3 ) and quadrangular 
flow (rq) are smaller than for elliptic flow (7-2), and there¬ 
fore harder to measure, they seem to be slightly better 
described by hydrodynamics with MC-KLN initial condi¬ 


tions than for MC-Glauber initial profiles. While the sen¬ 
sitivity to 77 /s is small, 7 + and 7-4 appear to give a slight 
preference to larger 77 /s values, consistent with the quali¬ 
tative conclusion drawn in the preceding section from the 
overall trend of the v n power spectrum. 

In summary, flow factorization breaking effects appear 
to open a window on the initial fluctuation spectrum that 
is only slightly blurred by uncertainties in the shear vis¬ 
cosity during the dynamical evolution from initial to fi¬ 
nal state. This observable therefore complements the 74- 
dependence of the magnitudes of the p T -integrated flow 
coefficients v n {2} for which shear viscosity effects domi¬ 
nate over variations in the initial-state fluctuation spec¬ 
trum. While neither the MC-Glauber nor the MC-KLN 
initial fluctuation spectrum allows to quantitatively re¬ 
produce all available data simultaneously, this observa¬ 
tion suggests that, by using the full complement of ex¬ 
perimentally accessible flow and fluctuation observables, 
it will in the future be possible to identify the correct 
initial-state model and, at the same time, quantify the 
quark-gluon plasma transport properties. 


V. CONCLUSIONS 

In this paper, we presented anisotropic flow studies 
for 0-0.2% ultra-central Pb+Pb collisions at the LHC, 
using the MC-Glauber and the MC-KLN initial condition 
models and evolving the system hydrodynamically with 
77 /s = 0, 0.08, 0.12, and 0.20. In the MC-Glauber model, 
we implement multiplicity fluctuations in the initial state, 
which boost the initial eccentricities e n {2} for n <9 by 
15-45%. A comparison with CMS data reveals that both 
MC-Glauber and MC-KLN models fail to reproduce the 
p-r-integrated v n hierarchy, especially the V 2 { 2 } /443 { 2 } 
ratio. Further comparisons with the pr-differential v n 
tend to favor a relatively large average value of 77 /s > 0.2 
for the medium. 

We found a large breaking of flow factorization for the 
elliptic flow coefficient in ultra-central collisions. Both 
the fluctuations of the flow magnitude and of the flow 
angle are important contributors to this breaking. Al¬ 
though our simulations can not fully reproduce the mea¬ 
sured v n spectrum, calculations of the r n agree overall 
quite well with the CMS measurements. All qualitative 
features and trends of the CMS data are correctly re¬ 
produced by the hydrodynamic model. Consistent with 
the conclusion from the v n comparison, the CMS r n data 
again slightly favor 77 /s ~ 0.2 over smaller values of the 
average specific shear viscosity during the dynamical evo¬ 
lution. 
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